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In this work we study the warm equation of state of asymmetric nuclear matter in the quark 
meson coupling model which incorporates explicitly quark degrees of freedom, with quarks coupled 
to scalar, vector and isovector mesons. Mechanical and chemical instabilities are discussed as a 
function of density and isospin asymmetry. The binodal section, essential in the study of the liquid- 
gas phase transition is also constructed and discussed. The main results for the equation of state 
are compared with two common parametrizations used in the non-linear Walecka model and the 
differences are outlined. 



I. INTRODUCTION 

The understanding of the properties of hot dense 
nuclear matter is an important problem in theoretical 
physics in the context of neutron stars as well as in 
high-energy heavy-ion collisions experiments where nu- 
clei undergo violent collisions [2(. Such collisions may 
produce states of nuclear matter at conditions far from 
those normally encountered in low energy collisions. For 
example, at high density or high temperature nucleons 
in the nucleus may dissolve into quarks and gluons. At 
low temperature, which can be attained in medium en- 
ergy heavy-ion reactions, no such exotic state can be pro- 
duced but there is the possibility of interesting liquid-gas 
phase transitions leading to the breakup of heated nu- 
clei into small clusters or droplets of nucleons Q. Such 
phase transitions have been identified in multifragmenta- 
tion experiments and in the crusts of neutron stars 0, Q . 
The possibilities of such phase transitions were previ- 
ously considered by several authors using different ap- 
proaches H0i|. 

In the present paper, we employ the quark-meson cou- 
pling model (QMC) [3, to investigate the liquid-gas 
phase transition in nuclear matter at different isospin 
asymmetries. In the QMC model, nuclear matter is 
described as a system of non-overlapping MIT bags 
which interact through the effective scalar and vector 
mean fields, very much in the same way as in the 
Walecka model, or more generally Quantum Hadrody- 
namics (QHD) Many applications and extensions of 
the model have been made in the last years [T3. IT3I Phil 
Il5| , including droplet formation at zero temperature |16| . 
An special improvement is related with a density depen- 
dent bag parameter. This model is known as modified 
QMC model (MQMC) [H G3 and some of its charac- 
teristics are discussed in the last section of the present 
paper. 

The crucial difference is that in the QMC, the effec- 
tive meson fields couple directly to the confined quarks 
inside the nucleon bags instead to point like nucleons as 
in QHD. While the QMC model shares many similari- 
ties with QHD-type models, it however offers new op- 



portunities for studying nuclear matter properties. One 
of the most attractive aspects of the model is that dif- 
ferent phases of hadronic matter, from very low to very 
high baryon densities and temperatures, can be described 
within the same underlying model. In the QMC, matter 
at low densities and temperatures is a system of nucleons 
interacting through meson fields, with quarks and glu- 
ons confined within MIT bags. For matter at very high 
density and/or temperature, one expects that baryons 
and mesons dissolve and the entire system of quarks and 
gluons becomes confined within a single, big MIT bag. 
Another important aspect of the QMC is that the in- 
ternal structure of the nucleon is introduced explicitly. 
Although the internal structure of nucleons can be effec- 
tively taken into account with an effective field theory 
employing hadronic degrees of freedom in a derivative 
expansion [19(, it is clear that such an approach will fail 
to describe the transition to a quark-gluon phase. It is 
found that the equation of state (EOS) for infinite nuclear 
matter at zero temperature derived from the QMC model 
is much softer than the one obtained in the Walecka 
model [llj. Also, the QMC nucleon effective mass lies in 
the range 0.7 to 0.8 of the free nucleon mass, which agrees 
with results derived from nonrelativistic analysis of scat- 
tering of neutrons from lead nuclei ppj | and is larger to 
the effective nucleon mass in the Walecka model. At finite 
temperature, there arises yet another difference between 
predictions of QMC and QHD, namely the behavior of 
the effective nucleon mass with the temperature at fixed 
baryon density. While in QHD-type models the nucleon 
mass always decreases with temperature, in the QMC it 
increases. The difference arises because of the explicit 
treatment of the internal structure of the nucleon in the 
QMC. When the bag is heated up, quark-antiquark pairs 
are excited in the interior of the bag, increasing the in- 
ternal energy of the bag. 

Here, we analyze the QMC model at finite tempera- 
ture for different proton fractions. The proper treatment 
of the problem with isospin asymmetry requires not only 
isoscalar scalar and vector mesons fields, but also isovec- 
tor meson fields. In particular we investigate the conse- 
quences of the different thermal behavior of the nucleon 
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mass in the QMC model on the liquid-gas phase transi- 
tion. The results are also compared with the ones ob- 
tained with two commonly used parametrizations for a 
non- linear Walecka model (NLWM), namely, NL-3 |2l| 
and TM-1 We organize the present paper as fol- 

lows. In Section II we describe the QMC model at finite 
temperature for asymmetric nuclear matter and discuss 
the binodal section. In Section III we present our numer- 
ical results and discuss differences with other approaches. 
Our summary is presented in Section IV. 



At finite temperatures, the three quarks inside the bag 
can be thermally excited to higher states and also quark- 
antiquark pairs can be created. For simplicity, we assume 
that the bag describing the nuclcon continues to remain 
in a spherical shape with radius R, which is now temper- 
ature dependent. The single-particle energies in units of 
R^ 1 are given as 



e™ = a?< + R(v u ±-v p ), 



(5) 



for the quarks and 



II. QMC AT FINITE TEMPERATURES 
A. Asymmetric matter at finite temperature 

Recently the QMC model has been generalized to fi- 
nite temperature taking a medium dependent bag con- 
stant |l-l|. In this generalized model, symmetric nuclear 
matter was considered and hence only the interaction of 
quarks through the exchange of effective isoscalar scalar 
(a) and vector (u>) mesonic fields were taken into account. 
We now extend this model to asymmetric nuclear matter 
at finite temperatures and include the contribution of the 
isovector vector meson p, in addition to a and u. 

In the QMC model, the nucleon in nuclear matter is 
assumed to be described by a static MIT bag in which 
quarks couple to effective meson fields, which are treated 
as classical in a mean field approximation. The quark 
field ip g (r, t) inside the bag then satisfies the Dirac equa- 
tion 



ij-d- (ro° - V a ) - 7°(K + K 3q V p ) 



^(r,t)=0, (1) 



where V a = g q a a , V u = g q Lo and V p = g q b 03 with er , 
loq and &03 being the classical meson fields. g q , g q and 
g q p are the quark couplings with the a, u> and p mesons 
respectively, mj is the current quark mass and r 3q is the 
third component of the Pauli matrices. In the present 
paper we consider only nonstrange q — u and d quarks 
only. 

At zero temperature, the normalized ground state (s- 
state) wave function for a quark in the bag is given as 



V> q (r, t) = Nexp(- 
where 

with 

and Xq is a Pauli spinor. 
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m = m - g q a a, 



(2) 



(3) 



(4) 



e - = - R<y u ± -v f 



(6) 



for the antiquarks, where the + sign is for u quarks and 
— for d quarks, and 



(7) 



The eigenvalues x nK for the state characterized by n and 
k are determined by the boundary condition at the bag 
surface, 



z 7 • m/C = C K - 



(8) 



Thus, the quark eigenvalues x nK become modified by the 
surrounding nucleon medium at finite temperature. 

The total energy from the quarks and antiquarks at 
finite temperature is 



Etot = ^2 (/I + /n«) > 



where 



and 



f q 



1 



f 

■I ! l . 



uk e (Q™/R+u q )/T + y 



(9) 



(10) 



(11) 



with v q being the effective quark chemical potential, re- 
lated to the true quark chemical potential fi q as 



Vq = fX q - V u - m q T V p . 

The bag energy now becomes 

Ebag = E to t — + -ttR 3 B, 
H o 



(12) 



(13) 



where B is the bag constant and Z parametrizes the sum 
of the center-of-mass (cm.) motion and gluonic correc- 
tions. Note that this center of mass correction is different 
from that of ^} • The effective nucleon mass is obtained 
from the bag energy and reads: 



MZr = Ebaq 



(14) 
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The bag radius R is determined by the equilibrium con- 
dition for the nucleon bag in the medium 



dM 



N 



OR 



= 0. 



(15) 



Once the bag radius is obtained, the effective nucleon 
mass is immediately determined. For a given tempera- 
ture T and scalar field a, the effective quark chemical 
potentials v q , q — u, d are determined from the total 
number of quarks of each type in the proton and in the 
neutron, i.e., 

4 = E (A - fU) s 3 > ( 16 ) 

q,n,K 

4 = E 2 "M«) (fn, - fl q ) = 2 ™r(j), ( 17 ) 
q,n,K 

for j = p, n. 

From the last two expressions one can see that the 
masses for protons and neutrons are different, since the 
effective chemical potentials v qi q = u,d are not the 
same. Nevertheless, it is a reasonable approximation, 
for the purposes of the present paper, to take for the nu- 
cleon mass the average value, so that we are left only with 
Eq. Ijltil) and hence, v u — v^. As a starting point, this 
approximation would go along the Walecka model which, 
with the most common parametrizations, does not distin- 
guish the masses of protons and neutrons. The situation 
would be different for problems where the neutron-proton 
mass difference is the main issue, such as those discussed 
in Ref . 23] . We stress that for a fixed T and a the quark 
distribution functions only depend on the effective chemi- 
cal potentials v q , which are determined by the constraints 
in Eqs. ((TBj) and (|TT|) . Therefore, just as it happens at 
zero temperature, also at finite temperature the effective 
mass does not depend on u> and &03- 

Before we proceed, we note that we are not considering 
density/temperature dependence in the zero point energy 
parametrized by Z because this would introduce extra, 
presently unkno wn p arameters. Blunden and Miller (see 
first reference in |13| ) have proposed changing Z through 
a linear dependence with the scalar field and concluded 
that for reasonable parameter ranges this has little effect. 

The total energy density at finite temperature T and 
at finite baryon density is 

£ = T^p E / d3k [ £ * (fi + fi) + v «(/< - /i)] 



± _2 2 1 2 2 1 2t.2 

—vn„o m, ,uj m n o n o, 

2 ^ 2 



(18) 



where /, and fi are the thermal distribution functions for 
the baryons and antibaryons, 



fi = 



1 



,(e*-ut)/T + I 



and ft = 



1 



; (e*+v ( )/T + I' 



(19) 



with i = p,n and e* = (k 2 + M* N 2 ) l/2 the effective nu- 
cleon energy, vi = /i; — Voi the effective baryon chemical 



potential and Voi = g^u + m T <7 P &03 ( m r = ±1/2 respec- 
tively for protons and neutrons). 

The thermodynamic grand potential density is defined 

as 



n 



£ -TS - E ViPi 



(20) 



with the entropy density S = S/V given by 
2 



S 



E 



(2tt) 3 



t—p,n 

+ /< In £ + (1-/0111(1-/0 



/ i ln/ i + (l-/0ln(l-/0 
(21) 



The proton or neutron density is given by 
2 



Pi 



(2tt) 3 



d 6 k (fi-fi), 



(22) 



so that the baryon density is p = p p + p n and the (third 
component of) isospin density p3 = p v — p n . The proton 
fraction is defined as 



Vp 



Pp 
P ' 



(23) 



The pressure is the negative of f2, which after an inte- 
gration by parts can be written as 



P 



3 ^ (2tt) 3 

i—p^n 



2 1 1 2 2 
—m„o~ H TO, ,Ll> 

2 2 u 



d 6 k- 



e*(k) 
1 



(fi + fi) 



2i2 

2 m P b 03- 



(24) 



From the above expression the pressure depends explic- 
itly on the meson mean fields a, u> and 6o3- It also de- 
pends on the nucleon effective mass Mjy which in turn 
also depends on the sigma field (see Eqs. 191141) ). 

At a given temperature and for given p p , p n , the effec- 
tive nucleon mass is known for given values of the me- 
son fields, once the bag radius R and the effective quark 
chemical potentials v q are calculated by using Eqs. ifFST) 
and (|16fl respectively. Maximizing the pressure with re- 
spect to the fields for a given temperature T and given 
chemical potentials pi , we obtain the following equations 
for the lu and 603 fields 



9u 



b()3 = 



9p 
2m 2 



A3- 



The <7 meson field is determined through 



dP 
&7 



dM 



dP 



dM* N J^ T da 



N 



dP 
~da~ 



(25) 



(26) 



MX 



We next calculate the binodal sections within this model. 
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B. Binodal section 

Nuclear matter is not stable at all temperatures, pro- 
ton fractions and densities. If the Gibbs energy of a two 
component phase is lower than the Gibbs energy of a 
one component phase, the system will separate into two 
phases [g, • The stability criteria may be expressed by 
the following relations 



du 
df 

dfip 



> 0, K 



p,v P 



dP 
dp 



> 0, (27) 



T,y p 



> and 



T,P 



dp n 

dy P 



<o, 



(28) 



T,P 



where u = E/B is the energy per particle. These condi- 
tions guarantee thermodynamical stability (positive spe- 
cific heat C v ), mechanical stability (positive compression 
modulus K) and diffusive stability (in a stable system, 
energy is required to change the proton concentration 
while the pressure and the temperature are kept con- 
stant). If any of these three criteria is violated there will 
occur a phase separation. 

The surface of two phase coexistence, the binodal 
surface, is determined imposing the Gibbs conditions, 
namely, for a given temperature, the pressure and chem- 
ical potentials are equal for both phases. A discussion 
of the properties of phase separation in different mod- 
els in terms of temperature, pressure and proton fraction 
is most conveniently done by comparing the respective 
binodal surfaces. Understanding the properties of the 
mixed phase of nuclear matter is important for different 
systems both in astrophysics and nuclear physics. One 
example is the crust composition of a neutron star, which 
near the surface is thought to form a Coulomb lattice of 
nuclei immersed in an electron gas and to be responsi- 
ble for the glitch phenomenon Also, in nuclear 
physics, highly excited nuclei created in heavy- ion col- 
lisions gives rise to multifragmentation which is inter- 
preted as a liquid-gas phase transition 

In order to obtain the proton and neutron chemical 
potentials in the two coexisting phases for a fixed pres- 
sure, we have used the geometrical construction [g, Q 
with the neutron and proton chemical potential isobars 
as a function of proton fraction. This takes into account 
the diffusive stability conditions Eq. (|28|l and the Gibbs 
conditions for phase separation. For a given temperature, 
the binodal section, which contains points under the same 
pressure for different proton fractions, is obtained from 
the above conditions, simultaneously with the following 
equations: 



P{v p ,v n ,M*) = P{v' p ,v' n ,M*'), 
Hi(i> p , v n , M*) = /j,i(v' p) v' n , M*'), i = p,; 



dP 



= 0, 



and 



dP 



= 0. 



(29) 



(30) 



(31) 



At this point, it is worth mentioning that the construc- 
tion of the binodal section for the QMC is slightly more 
complicated than for the non-linear Walecka model. This 
is because Eq. l|15fl has to be enforced in the numerical 
calculation for every a. 

The binodal section is formed by two branches, one cor- 
responding to a gas phase and small proton fraction and 
the other to a liquid phase and large proton fraction. It 
is, therefore, energetically favorable for nuclear matter to 
separate into a liquid phase with a large proton fraction 
(less asymmetric) and a gas phase with a small proton 
fraction (more asymmetric). This behavior is common 
to other relativistic mean-field models and reflects the 
fact that the contribution from the p-meson gives a term 
in the energy per particle of the form e sym (pB)8 2 where 
the asymmetry parameter is <5 — (p n — p p )/ pb and the 
coefficient e sym increases with ps- 

In the next section we will display the binodal section 
results and compare the same with different non-linear 
relativistic models. 



III. RESULTS 

We start by fixing the free-space bag properties. We 
have used zero quark masses only and Rq = 0.6 fm for 
the bag radius. There are two unknowns, Z and the bag 
constant, B. These are obtained as usual by fitting the 
nucleon mass, M = 939 MeV and enforcing the stability 
condition for the bag. The values obtained for Z and 
B are displayed in Table 1. The quark-meson coupling 
constants g%, g^ — 3g* and g p — g q p are fitted to obtain 
the correct saturation properties of nuclear matter, Eb = 
e/p - M = -15.7 MeV at p = p Q = 0.15 fm" 3 and 
o-sym = 32.5 MeV. We take the standard values for the 
meson masses, also shown in Table 1. 

In the sequel we will frequently compare QMC to 
two different parametrizations of the non-linear Walecka 
model (NLWM), namely the NL-3 parametrization with 
non-linear contributions from the scalar meson only [2l| 
and the TM-1 parametrization which includes quartic 
terms on scalar and vector meson |22(. 

We first plot the energy per baryon as a function of the 
baryon density in Fig. ^ for symmetric nuclear matter, 
y p = 0.5, for temperatures 0, 4, 8, 12, 16 and 20 MeV. 
As expected, this functional at zero temperature has a 
minimum at the nuclear matter saturation density, po, 
corresponding to a binding energy of —15.7 MeV. As the 
temperature increases the minimum shifts towards higher 
densities. This may be understood from the fact that to 
compensate for the larger kinetic energy a larger value 
of p is needed to give the minimum. For larger densities 
the nuclear repulsion effects take over and again energy 
increases. For higher temperatures, the minimum of the 
curves become positive, as in the usual NLWM. In Fig. 21 
we plot the binding energy for neutron matter for tem- 
peratures 20, 50 and 100 MeV and compare the results 
with TM-1 and NL-3 parametrizations of the NLWM. 
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At higher temperatures, the QMC model gives a larger 
binding energy as compared to the Walecka models. The 
reason for this is the contribution for the thermal energy 
of the nucleon bag which is absent in the Walecka model. 

This can be clearly seen from Fig. where the baryon 
effective masses are plotted as a function of density for 
different temperatures. In the same figure, we have also 
shown the results of the parameter set NL-3 and TM-1 
for comparison. The value of M* increases with temper- 
ature. We have checked that at zero density and tem- 
perature ~ 40 MeV, the QMC effective mass increases 
from 939 MeV in contrast to the calculation of hot nu- 
clear matter using Walecka model. The reason for this is 
that the sigma field is not strong at higher temperatures. 
Also, there is a significant contribution to the effective 
mass arising from the thermal excitation of the quarks 
inside the nucleon bag, which adds to the mass of the 
nucleon at higher temperatures. This contribution which 
is absent in the Walecka model calculation, appears dom- 
inant over the contribution from the sigma field and the 
net effect is a rise of the effective mass. The behavior of 
the effective mass at high temperature obtained here is 
contrary to the results presented in Ref. 0] , where tem- 
perature was introduced only at the hadron level, and 
therefore the behavior of the effective mass with tem- 
perature not-surprisingly is equivalent to the results of 
Walecka-type of models. 

In 2] we plot the effective radius in units of the free nu- 
cleon radius i? = 0.6 fm as a function of the baryon den- 
sity for temperatures 20, 50 and 100 MeV. It is obesrved 
that the nucleon bag shrinks with the increase of tem- 
perature. A similar behaviour is also predicted with the 
MQMC model [3. Such an effect is also encountered 
using the thermal skyrmion |25| where a baryon shrinks 
with temperature. We have also calculated the effective 
mass of the nucleon as a function of density at differ- 
ent temperatures for the QMC and the MQMC models. 
As already discussed in 0, ^ the effective mass in the 
MQMC decreases faster with p than in the QMC model. 
However, also in the MQMC there is and increase of the 
effective mass with temperature. A comparison between 
the QMC model and different versions of the modified 
MQMC has been performed in Ref. • O ne important 
conclusion was that, contrary to the QMC result, the 
bag radius increases with density for all MQMC mod- 
els. For densities not much larger than nuclear matter 
saturation density po the bags start to overlap which im- 
plies a breakdown of the model. On the other hand the 
MQMC models contain sufficient free parameters which 
allow a good reproduction of the ground state properties 
of nuclei. In the present work we were interested in de- 
scribing propeties of nuclear matter for densities which go 
beyond the saturation density and therefore have chosen 
to consider the QMC model. 

The possible existence of a liquid-gas phase transition 
is determined by the pressure. We plot this quantity 
as a function of the baryon density, p, for low temper- 
atures in Fig. [S] for symmetric nuclear matter. At zero 



temperature, the pressure decreases with density, reaches 
a minimum, then increases and passes through P = 
at p = po, where the binding energy per nucleon is a 
minimum. A decrease of the pressure with density corre- 
sponds to a negative compressibility, K = 9(dP/dp), and 
is a sign of mechanical instability. As the temperature in- 
creases the region of mechanical instability decreases. At 
T=17.7 MeV, the pocket disappears. This corresponds 
to the critical temperature defined by {dP/dp)T, Vp = = 
(d 2 P/dp 2 )T. Vp = 0, above which the liquid-gas phase 
transition is continuous. It is comparable to the values 
for the critical temperature obtained with Skyrme inter- 
actions [2(j or relativistic models 0. 

In Fig. El we plot the pressure as function of density 
for temperatures 20, 50 and 100 MeV. In the same fig- 
ure we have also shown the results obtained with the 
NLWM of TM-1 and NL-3 parameter sets to compare. It 
is found that the QMC model gives a softer EOS, which is 
true also for any proton fraction. Although QMC gives a 
higher energy per particle than NL-3 and TM-1, the pres- 
sure, which corresponds to the derivative of the energy 
per particle with respect to density, does not increase so 
fast with p. 

In Fig. [7| the pressure is given as a function of the 
baryon density, p, for different temperatures and differ- 
ent proton fractions. The region of mechanical instability 
decreases both with the increase of temperature and the 
decrease of the proton fraction. This is clearly seen in 
Fig. [SI where the quantity (dP/dp)T. Vp is shown for two 
different temperatures and three different proton frac- 
tions. The curve corresponding to T = 10 MeV and 
y p = 0.3 has the smallest region with a negative value of 
the derivative (dP / dp)x,y v - 

The behavior of the critical temperature with the pro- 
ton fraction is shown in Fig. It decreases rapidly with 
y p and is going to zero for y p = 0.041. Although there 
will be no mechanical instabilities above T c , chemical in- 
stabilities can still develop above T c . The present results 
are comparable to the ones discussed in Refs. 0E3jI13- 

The entropy per baryon is an important quantity to 
describe a collapsing massive star. During the whole pro- 
cess of the collapse the entropy per baryon, including the 
contribution of baryons and leptons, remains low: it is 
of the order of S/p ~ 1 initially, increases slightly before 
neutrino trapping and remains constant, S/p ~ 1 — 2, 
afterwards [29]. According to different models the con- 
tribution from the baryons alone changes from 80-90% at 
the saturation density to 50-65% at 10 times that density 
0. In Fig. llOl the entropy per baryon is plotted and com- 
pared with the corresponding quantity calculated within 
the TM-1 parametrization of NLWM. We conclude that 
QMC gives values compatible with the ones discussed 
in Ref. pj. The entropy per baryon in QMC decreases 
slower than in TM-1. 

In Fig. 11, we plot the scalar and vector potentials as 
a function of density for T=8 MeV and y p — 0.4. In the 
same figure we have also shown the non-linear Walecka 
model results. The ui and a contributions in the QMC 
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are both weaker than the corresponding contributions in 
the NL-3 and TM-1 parametrizations of the NLWM. As 
a consequence, the QMC effective mass changes less with 
density than the TM-1 and NL-3 effective masses. This 
can also be confirmed in Fig. 3. A second consequence is 
a softer EOS due to the weaker omega field. 

Isobars of nuclear pressure are plotted on the p, p and 
fi n space in Fig. 12 for T = 6 MeV for two different 
pressures, P = 0.1 and 0.063 MeV/fm 3 . The diffusively 
unstable regions can be seen clearly in chemical poten- 
tial isobars in this figure. According to the inequality 
|[S8J, the region of negative (positive) slope for p p (fi n ) 
is unstable. Violation of the stability criteria is an indi- 
cation of phase separation. The surface of the two phase 
coexistence in the {p, T, y p ) space, the binodal surface, 
is obtained from conditions H29fl - I|31|) . I n Fig. 13 we 
plot the binodal sections for the QMC and the NL-3 and 
TM-1 parametrizations, at T= 6 MeV. The critical point 
occurs for similar values of the pressure and proton frac- 
tion in all three models. For the TM-1 parametrization, 
the critical pressure is 10% larger than the QMC criti- 
cal pressure. On the other hand the proton fraction for 
NL-3 is ~ 10% larger than the corresponding values for 
QMC. The largest difference between the models is in the 
shape: the region of configuration space where the phase 
separation is favorable is smaller in the QMC and occurs 
at smaller values of y p in the liquid phase. As a con- 
sequence chemical instability occurs for smaller isospin 
asymmetries, which may have consequences in isospin 
fractionation |28| . The finite temperature results agree 
with the ones obtained within QMC at zero temperature, 
16]. For larger temperatures, the binodal has a similar 
shape but extends to smaller isospin asymmetries Q, || . 

IV. SUMMARY 

We now summarize the results and conclusions of the 
present work. We have studied asymmetric nuclear mat- 
ter at finite temperature using the QMC model which 
incorporates explicitly quark degrees of freedom. The re- 
sults for the EOS were compared with the two non-linear 
Walecka models, namely NL-3 and TM-1 parametriza- 
tions. The mean effective fields a, lu and p are determined 
from the minimization of the thermodynamical potential, 
and the temperature dependent effective bag radius was 



calculated from the minimization of the effective mass of 
the nucleon mass of the bag. The thermal contributions 
of the quarks, which are absent in the Walecka model, 
arc dominant and lead to a rise of the effective nucleon 
mass at finite temperature. This is contrary to the results 
presented in Ref. [l5( , where temperature was introduced 
only at the hadron level, and therefore the behavior of 
the effective mass with temperature is equivalent to the 
results of Walecka- type models 11]. In the present cal- 
culation, the effective radius of the nucleon bag shrinks 
with increase of the temperature. 

The equation of state as derived here is softer than 
the ones obtained within the non-linear Walecka model 
for the NL-3 and TM-1 parametrizations, for all temper- 
atures. The region of mechanical instability decreases 
with the increase of T and decrease of y p . Also, the criti- 
cal temperature decreases with the proton fraction, from 
a maximum value of T = 17.7 MeV at y p = 0.5 to T = 
MeV at y p = 0.041. We have shown that the potentials in 
the QMC are weaker than the corresponding ones in the 
NLWM for the NL-3 and TM-1 parametrizations, which 
imply a softer EOS and a weaker change of the effective 
mass with p. 

We have also constructed and studied the binodal sec- 
tions which is essential for studying the liquid-gas phase 
transition and to understand under which conditions 
chemical instabilities develop, leading to isospin fraction- 
ation. It is clear from the binodal that the system prefers 
to separate into regions of higher density and smaller 
isospin asymmetry and regions of lower density and large 
isospin asymmetry 0,H,Hg. The phase separation in the 
QMC occurs at smaller isospin asymmetry, correspond- 
ing to lower values of y p , than the ones predicted by the 
NL-3 and TM-1 parametrizations of the NLWM. 

Extensions of the formalism presented here to droplet 
formation are currently under process. 
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TABLE I: Parameters used in the calculation for Rq = 0.6 
fm. All masses and bag pressure are in MeV. 



m q B 1 ' 4 


Z 


m CT 


m w 


m p 


3u> 


gl 


9p 


211.3 


3.987 


550 


783 


763 


8.9539 


5.981 


8.655 




P(fm 3 ) 



FIG. 1: Binding energy per nucleon f?jasa function of nu- 
clear matter density, p, at temperatures T = 0, 4, 8, 12, 16 
and 20 MeV. 




FIG. 2: Binding energy per nucleon Eb as a function of nu- 
clear matter density, p, for neutron matter at temperatures T 
= 20, 50 and 100 MeV respectively. 




FIG. 3: Effective mass of the nucleon versus density at tem- 
peratures 20, 50 and 100 MeV in QMC, NL-3 and TM-1 pa- 
rameters. 
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FIG. 4: The effective bag radius as a function of nuclear mat- 
ter density p for temperatures T — 20, 50 and 100 MeV re- 
spectively. 




FIG. 5: Pressure P as a function of nuclear matter density, 
p, for symmetric matter at temperatures T = 0, 4, 8, 12, 16 
and 20 MeV respectively. 




FIG. 6: Pressure P as a function of nuclear matter density, p, 
for neutron matter at temperatures T = 20, 50 and 100 MeV, 
respectively, from bottom to top. Note that the QMC model 
gives softer equation of state. 




FIG. 7: Pressure P as a function of nuclear matter density, p, 
at different temperatures and different proton fractions. 
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FIG. 8: dP/dp as a function of the proton fraction. 




FIG. 9: Critical temperature as a function of the proton frac- 
tion. 




FIG. 10: Entropy per baryon, S/B, calculated within QMC 
and TM-1, as a function of nuclear matter density, p, for 
symmetric matter at temperatures T = 10, 20 and 40 MeV 
from bottom to top. 




FIG. 11: Potentials versus nuclear matter density, p, at tem- 
perature = 8 MeV and y p = 0.4 in QMC, NL-3 and TM-1 
parameters. 



13 




FIG. 12: Nuclear pressure isobars for the proton (lower) and 
neutron (upper) chemical potentials as a function of proton 
fraction for P = 0.1 and 0.063 MeV/fm 3 . 




FIG. 13: The binodal section at T=6 MeV for QMC, NL-3 
and TM-1 



